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ABSTRACT 

We investigate the electrodynamic structure of a pulsar outer-magnetospheric particle accelerator 
and the resultant gamma-ray emission. By considering the condition for the accelerator to be self- 
sustained, we derive how the trans-magnetic-field thickness of the accelerator evolves with the pulsar 
age. It is found that the thickness is small but increases steadily if the neutron-star envelope is 
contaminated by sufficient light elements. For such a light element envelope, the gamma-ray luminosity 
of the accelerator is kept approximately constant as a function of age in the initial ten thousand years, 
forming the lower bound of the observed distribution of the gamma-ray luminosity of rotation-powered 
pulsars. If the envelope consists of only heavy elements, on the other hand, the thickness is greater 
but increases less rapidly than what a light element envelope has. For such a heavy element envelope, 
the gamma-ray luminosity decreases relatively rapidly, forming the upper bound of the observed 
distribution. The gamma-ray luminosity of a general pulsar resides between these two extreme cases, 
reflecting the envelope composition and the magnetic inclination angle with respect to the rotation 
axis. The cutoff energy of the primary curvature emission is regulated below several GeV even for 
young pulsars, because the gap thickness, and hence the acceleration electric field is suppressed by 
the polarization of the produced pairs. 

Subject headings: gamma rays: stars — magnetic fields — methods: analytical — methods: numerical 
— stars: neutron 



1. INTRODUCTION 

The Large Ar ea Telescope aboard Fermi Gamma-ray 
Space Telescope (jAtwood et al. 2009) has proved remark- 
ably successful at discovering rotation-powered pulsars 
emitting photons above 0.1 GeV. Thanks to its superb 
sensitivity, the number of gamma-ray pulsars has in- 
creased f rom six in Compton Gamma Ray Observa- 
tory era ([Thompson! |2004) to more than one hundred 
(Nolan 2012|). Plotting their best estimate of the gamma- 
ray luminosity, L 7 , against the spin-down luminosity, 
L spin = iir 2 IPP~ 3 , lAbdo et al.l pOlOh found the im- 
portant relation that L 7 is approximately proportional 
to i sp in ' 5 (with a large scatter), where I refers to the 
neutron-star (NS) moment of inertia, P the NS rota- 
tional period, and P its temporal derivative. How- 
ever, it is unclear why this relationship arises, in spite 
of its potential importance to discrimina te pulsar emis- 
sion models such as the pola r -cap model (Hardin g et al.l 
119781 : iDaughertv fc Harding! [l982t iDermerl I1994T) . the 
outer-gap model (I Chiang fc Romanilll992HRomanij 
Zhang fe Chengill997t iTakata et al. 2004t lfflroteni 



1996 



2008 



Wang et aLl 120111) the pair-st arved polar-cap mode l 
(|Venter et all 12009( 1 (see also lYuki fc Shibatal (poll ) 
for the possible co-existence of such models), and 



the emission model from the wind zone (IPetrl 12011 
iBai fc Spitkovskll2010al lb1; lAharonian et al.ll2012l l 

Recent gamma-ray observations suggest that the 
pulsed gamma-rays are emitted from the higher alti- 
tudes of a pulsar magnetosphere. This is because the 

1 Postal address: TIARA, Department of Physics, National Ts- 
ing Hua University, 101, Sec. 2, Kuang Fu Rd.,Hsinchu, Taiwan 
300 



observed light curves (|Abdo et al.l [2010) favor fan-like 
emission geometry, which scan over a large fraction of 
the celestial sphere, and because the Crab p ulsar shows 
pulsed photons near and above 100 GeV (|Aliu et al.l 
l2011tlAleksic et al.ll20liaHbT ). which rules out an emission 
from the lower altitudes, where strong magnetic absorp- 
tion takes place for 7-rays above 10 GeV. Consequently, 
higher- altitude emission m odels such as the outer-gap 
model ( Cheng et al.lll986alfbh, the h igh-altitude slot-gap 
model (|Muslimov fc Har ding 20(3, or the pair-starved 
polar-cap model ( Venter et a l. 2009) , gathered attention. 
It is noteworthy that the outer-gap model is presently the 
only higher-altitude emission model that is solvable f rom 
the basic equations self-consistently (Hirotani 201 la]). In 
the present paper, therefore, we focus on the outer-gap 
model and derive the observed relationship L~ ( oc L sp j n ' 5 
both analytically and numerically. 

We schematically depict the pulsar outer- 
magnetospheric accelerator (i.e., the outer gap) in 
figure m As the NS rotates, there appears the light 
cylinder, within which plasmas can co-rotate with the 
magnetosphere. The magnetic field lines that become 
tangential to the light cylinder at the light cylinder 
radius, tolc = cP/2ir, are called the last-open magnetic 
field lines, where c refers to the speed of light. Pairs 
are produced via photon-photon pair production mostly 
near the null-charge surface and quickly polarized by 
the magnetic-field-aligned electric field, Eh , in the 
gap. In this paper, we assume that the rotation and 
magnetic axes reside in the same hemisphere to obtain 
En > 0, which accelerates positrons (e + 's) outwards 
while electrons (e~'s) inwards. These ultra-relativistic 
particles have Lorentz factors, 7 
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Fig. 1. — Side view of an outer gap. The neutron star (filled 
circle on the left) obliquely rotates around the vertical axis with 
magnetic inclination angle a. The thin solid curves denote the 
magnetic field lines, while the dashed straight line the null-charge 
surface, on which the magnetic field lines become perpendicular 
to the rotation axis. Outside the light cylinder (the vertical long 
dashed line) , plasmas that are frozen to the magnetic field lines can 
only migrate outwards (as a pulsar wind) because of the causality 
requirement in special relativity. The light-cylinder radius, tolCj 
becomes typically a few or several hundred neutron-star radii for 
young pulsars. In modern outer-gap model s (Hirotani ct aLl l2003l ; 
Takata ct al. 2004; Hirotani 2006a; Ta kata et al.1120041) , it isproved 
that the outer gap extends between the stellar surface (because of 
a negative charge density in the lower altitudes) and the vicinity of 
the light cylinder (because of a positive charge density in the higher 
altitudes). Thus, typical inward (or outward) photons propagate 
distance Zi (or I2) before escaping from the gap (shaded region). 



photons efficiently by the curvature process. 

2. ANALYTICAL EXAMINATION OF OUTER-GAP 
LUMINOSITY 

In this section, we analytically derive how the gamma- 
ray luminosity of an outer gap evolves with time. In the 
outer magnetosphere, only the dipole component remains 
in the magnetic field; thus, the inhomogeneous part of 
the Maxwell equation (i.e., the Poisson equation for the 
electro-static potential) give s the magnetic-field-aligned 
electric field fHirotanill2008l ). 



E 



2n7 LC 



(1) 



where p denotes the NS magnetic dipole moment, and h m 
the trans-magnetic-field thickness of the gap. Since the 
Poisson equation is a second-order differential equation, 
E\\ is proportional to h m 2 . Electrons (e~'s) and positrons 
(e + 's) are created via photon-photon (and sometimes via 
magnetic) pair production, being subsequently polarized 
by and accelerated in the opposite directions, to fi- 
nally attain the terminal Lorentz factor 
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where p c refers to the radius of curvature of particle's 
motion in the three-dimensional magnetosphere, e the 
charge on the positron. Photons are radiated by such 
ultra-relativistic e s via curvature process with charac- 
teristic energy, 



hv r 
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(3) 



where h denotes the Planck constant, H = h/2%. Once 
h m is obtained, we can readily compute the 7-ray lu- 
minosity of cu rvature radiation from an outer gap by 
(|Hirotariill2008h 



L 1 a 2.36(i/F„) pcak x 4nd 2 f n « 1.23 fa^ 



p 2 n 4 



(4) 



where /n, which has been conventionally assumed to be 
approximately unity, refers to the flux correction factor 
(jRomani. R. fc Watterd 120101) , and fl = 2tt/P the rota- 
tion angular frequency of the NS. Here, it is assumed 
that the current density flowing i n the gap is compa- 
rable to the Goldreich-Julain value (jGoldreich. fc Julian! 
1969). The last factor, ^ 2 f2 4 /c 3 is proportional to the 
spin-down luminosity, £ S pm- Therefore, the evolution 
law, L 1 cc L®p in , is crucially governed by the evolution of 
h m as a function of the NS age, t. 

The evolution of h m is essentially controlled by the 
photon-photon pair production in the pulsar magneto- 
sphere. To analytically examine the pair production, we 
assume the static dipole magnetic field configuration for 
simplicity, and consider the plane on which both the ro- 
tational and magnetic axes reside (fig. [T]). On this two- 
dimensional latitudinal plane, the last-open field line in- 
tersects the NS surface at magnetic co-latitudinal angle 
0max ( measurec i from the magnetic dipole axis) that sat- 
isfies 



sin 



S11T 



(0LC-C*) 



otlc/ sin 6» LC 



(5) 



where r* denotes the NS radius, #lc the angle (measured 
from the rotation axis) of the point where the last-open 
field line becomes tangential to the light cylinder, and 
a the inclination angle of the dipole magnetic axis with 
respect to the rotation axis. A magnetic field line can be 
specified by the magnetic co-latitude (measured from the 
dipole axis) , , where it intersects the stellar surface. A 
magnetic field line does not close within the light cylinder 
(i.e., open to large distances) if < 6* < fl™ 8 *. Thus, the 
last-open field lines, 6* — #™ ax , corresponds to the lower 
boundary, which forms a surface in a three-dimensional 
magnetosphere, of the outer gap. 

Let us assume that the gap upper boundary coin- 
cides with the magnetic field lines that are specified by 
(?» = (1 — ft. m )0" lax - Numerical examinations show that 
h m , indeed, changes as a function of the distance along 
the field line and the magnetic azimuthal angle (mea- 
sured counter-clockwise around the dipole axis). Nev- 
ertheless, except for young pulsars like the Crab pulsar, 
an assumption of a spatially constant h m gives a rela- 
tively good estimate. Thus, for an analytical purpose, 
we adopt a constant /i m in this analytical examination. 
In this case, we can specify the middle-latitude field line 
by the magnetic co- latitude f?„ = (1 — ft, m /2)6'" lax . Screen- 
ing of En due to the polarization of the produced pairs, 
takes place mostly in the lower altitudes. It is, therefore, 
appropriate to evaluate the screening of E\\ around the 
point (ro,#o) where the null-charge surface intersects the 
middle-latitude field line (fig. [TJ . 

An inwardly migrating electron (or an outwardly mi- 
grating positron) emits photons inwards (or outwards), 
which propagate the typical distance l\ (or I2) before es- 
caping from the gap. Denoting the cross section of the 
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inward (or outward) horizontal line from the point (ro,#o) 
and the upper boundary as (ri,#i) (or as (r2,#2)), and 
noting cos 9$ = ri cos 9\ = T2 cos 02 , we obtain 



/ 1 = ro cos 6q (tan 9q — tan 0\ ) , 



(6) 



h = tq cos (tan 02 — tan 0q). (7) 
Along the upper-boundary field line, we obtain 
sin 2 (6»i - a) sin 2 (6» 2 - a) sin 2 [(1 - h m )df ax ] 



n 



(8) 



whereas along the middle-latitude field line, we obtain 



sin 



%-a) sin 2 [(l-V/2)0r x ] 



(9) 



Combining these two equations, and noting 0™ ax <g; l 7 
we find that 0\ (< 0q) and 82 (> #0) can be given by the 
solution that satisfies 



cos sin 



a) = 



1 - h m 

1 - 



cos 6*o sin 2 (^o — 01), 



where 9q is given by 
1 



tan #o=2'3 tan ol - 



\/9 tan 2 



(10) 



(11) 



Thus, if we specify a, we can solve 9 = 9\ and 9 — 9 2 as 
a function of /i m by equation pH|) . Substituting these #1 
and 6*2 into equations © and ([7]), we obtain l\ and fo, 
where ro depends on P. 

If /i m <C 1, we can expand the left-hand side of equa- 
tion (1101) around 9 — 9o, where 9 = 9\ for inward (or 
9 — 92 for outward) 7-rays to find 02 — #0 = #0 — 9\ oc 
\Jh m . That is, the leading terms in the expansion vanish 
and we obtain l\ = l 2 from the next-order terms, which 
are quadratic to 9 — 9q. Assuming Lx oc t~@, where 
(3 pa 0.48 is appropriate for t < 10 4 years for a light- 
element-envelope NS and for t < 10 5 years for a heavy- 
element NS, we find h m cx P 5 /V~ 1/6 ^ /2 , and hence 
L 7 oc p- 3 / 2 ^ 3 / 2 ^ 3 / 3 / 2 . Since the dipole radiation formula 
gives P oc ^i 1 / 2 , we obtain L 1 oc ^ 3 (/3-V 2 )/ 2 « £-0-0 3 . 
Thus, when the gap is very thin, which is expected for 
a light-element-envelope NS, L 7 little evolves with the 
pulsar age, t. 

However, if h m > 0.2, say, the rapidly expanding mag- 
netic flux tube gives asymmetric solution, I2 > h- That 
is, the third and higher order terms in the expansion con- 
tribute significantly compared to the quadratic terms. 
Thus, we must solve equation (fTU|) for 9 (= 9\ or 9 2 ) 
without assuming \9 — 0q\ <C 1, in general. 

Let us now consider the condition for a gap to be self- 
sustained. A single ingoing e~ or an outgoing e + emits 



(N, 



7J1 



or 



(TV. 



7J2 



eEuh/(hv c ) 



eE«l 2 /{hv c ) 



(12) 



(13) 



photons while running the typical distance l\ or I2, re- 
spectively. Such photons materialize as pairs with prob- 
ability 

n = ZiFicri/c (14) 



T2 = I2F2CF2/C, 



(15) 



where Pi and P2 denote the X-ray flux inside and outside 
of (ro,6>o), respectively; ui and o\ are the pair-production 
cross section for inward and outward photons, respec- 
tively. Thus, a single e" or e + cascades into 



pairs or into 



(iV 7 ) 2 T 2 



eP|| P 2 2 

< 2 0"2 



(16) 



(17) 



pairs within the gap. That is, a single inward-migrating 
e~ cascades into pairs with multiplicity (7V 7 )iTi. Such 
produced pairs are polarized by En. Each returning, 
outward-migrating e + cascades into pairs with multiplic- 
ity (iV 7 )2T2 in outer magnetosphere. As a result, a single 
inward e~ cascades eventually into (JV 7 )irx • (iV 7 ) 2 T 2 in- 
ward e _ 's, which should become unity for the gap to be 
self-sustained. Therefore, in a stationary gap, the gap 
thickness h m is automatically regulated so that the gap 
closure condition, 



(AT 7 ) 1 r 1 • (Ag 2 T 2 = 1, 



(18) 



may be satisfied. 

Approximately speaking, a single, inward-migrating e~ 
emits (-/V 7 )i ~ 10 curvature photons, a portion of which 
head-on collide the surface X-ray photons to material- 
ize as pairs with probability T\ ~ 10 -3 within the gap. 
Thus, a single e~ cascades into (iV 7 )iri ~ 10 pairs in 
the gap. Each produced e + return outwards to emit 
(Nj) 2 ~ 10 5 photons, which materialize as pairs with 
probability r 2 ~ 10~ 6 by tail-on colliding with the sur- 
face X-rays. In another word, (iV 7 )iTi ~ 10 holds regard- 
less of the nature of the pair production process (e.g., 
eithe r photon-photon or magnetic process (jTakata et al I 
2010)) in the lower altitudes, because it is determined 
by the pair-production efficiency in the outer magneto- 
sphere (iV 7 )2T2 ~ 0.1, which is always due to photon- 
photon pair production. 

In general, (iV 7 )i, t\, (N^)?, T2 are expressed in terms 
of h m , P, jjL, T, and a, where T denotes the NS surface 
temperature. Note that we can solve P = 2-k/Q. as a 
function of the NS age, t, from the spin-down law. Thus, 
specifying a and the cooling curve, T = T(t), we can 
solve h m as a function of t from the gap closure condition, 
(A^ 7 )iTi(7V 7 )2T2 = 1. Note also that the spin-down law 
readily gives the spin-down luminosity, L spin oc PP~ 3 , 
as a function of t, once P = P(t, a) is solved. On these 
grounds, L 7 can be related to L S pin with an intermediate 
parameter t, if we specify the cooling curve and the spin- 
down law. 

Substituting equations (fi"6j) and (fl"7|) into (fT8|) . we ob- 
tain 



ePii y/FiaiF 2 (T2 



1. 



(19) 



where 



F m = tt(1 - Mi) (-) f 



BAT) 
hv 



ap(v, v~!,Hi) (20) 
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with i — 1,2; v 1 denotes the 7-ray frequency, and B V (T) 
the Planck function. We have to integrate over the soft 
photon frequency v above the threshold energy 



h v. 



th,, 



2(m e c 



2 \2 



(1 - m)hv~, 



(21) 



where m e c 2 refers to the rest-mass energy of the electron. 
The cosine of the collision angle /Xj becomes 1 — /ii = 
1 — sin 8q for outward (or 1 — fj,2 = 1 + sin #0 for inward) 
7-rays. That is, collisions take place head-on (or tail-on) 
for inward (or outward) 7-rays. The total cross section 
is given by 



cr P = — er T (l - v 2 ) 
16 



( 3 _^ )ln i±^ .,.„., .,2., 
1 — v 



2v(2-v 2 ) 



where ctt denotes the Thomson cross section and 



(m e c 2 ) 2 



1 — ^ hvhv 1 



(22) 



(23) 



Pair production takes place when the 7-rays collide 
with the surface X-rays in the Wien regime, that is, at 
hp 3> kT. An accurate evaluation of 02 requires a careful 
treatment of the collision geometry, because the thresh- 
old energy, hvtb,2, strongly depends on the tiny collision 
angles. In the numerical method (next section), the pair- 
production absorption coefficient is explicitly computed 
at each point in the three-dimensional pulsar magneto- 
sphere by equation (|50[) . However, in this section, for 
analytical purpose, we simply adopt the empirical rela- 
tion, 

y/FiaiF 2 a 2 = ey/l — /zicrr-Fx, (24) 

where e = 0.004, 0.01, and 0.038 for a = 45°, 60°, and 
75°, respectively; 1 — fii « 2. The X-ray flux is evaluated 
at (ro^o) sucri that 



Fx = 



1 



2.82/cT 47rr 2 ' 



(25) 



where Lx refers to the luminosity of photon radiation 
from the the cooling NS surface. For a smaller a, the 
point (r^fii) is located in the higher altitudes, where the 
magnetic field lines begin to collimate along the rota- 
tion axis, deviating from the static dipole configuration. 
Thus, the collision angles near the light cylinder, and 
hence a 2 decreases with decreasing a. The explicit value 
of e can be computed only numerically, solving the pho- 
ton specific intensity from infrared to 7-ray energies in 
the three-dimensional pulsar magnetosphere. 

The last factor, l\l 2 i in the left-hand side of equa- 
tion (|19p is given by 

hh = ro 2 cos 2 8 (t&n9 - tan#i)(tan# 2 - tan# ) (26) 
Thus, equation (fH?t gives 

eE\\ L x /c r 

hv c 2.82/cT V ^ 

x cos 2 6»o (tan O - tan 9i ) (tan 9 2 - tan ) = 1, (27) 

where the tq dependence vanishes. Substituting equa- 
tions m, ©i © into (f27|) . we can solve h m as a function 
of Lx/kT, P, and /z. 



To describe the evolution of P 
adopt in this paper 



P(t) = 27r/fi(i), we 



2n4 



mn = c 



(28) 



where C = (2/3) sin 2 a for a magnetic dipole braking, 
while C = 1 + sin 2 a for a force- free braking ( Spitkovski 
2006). Assuming a magnetic dipole braking, we obtain 



P = 39.2ms/i 30 / 45 1/2 (t/10 3 years) 



!/ 2 , (29) 

where p, 30 = /it/(10 30 G cm 3 ) and J 45 = 7/(10 45 g cm 2 ). 
Thus, if we specify a cooling scenario, T = T(t), equa- 
tion (l27l) gives /i m as a function of t. Note that the a 
dependence of the spin-down law is not essential for the 
present purpose; thus, C — 2/3 is simply adopted. Once 
hm. — h m (t) is obtained, equation Q readily gives L 1 
as a function of i, and hence of L sv - ln . It is worth not- 
ing that the heated polar-cap emission is relatively weak 
compared to the cooling NS emission, except for millisec- 
ond or middle-aged pulsars. 

Let us now consider the cooling scenario. Since 
the mass of PSR J1614-2230 is precisely measured to 
be 1.97M Q (i.e., 1.97 solar masses), and since other 
three NSs (4U1700-377, B 1957+20, and J1748-2021B) 
(|Lattimer fc Prakashll200ll ) are supposed to be heavier 
than 2.0M Q , we exclude the equation of state (EOS) 
that gives smaller maximum mass than 1.95M0. That 
is, we do not consider the fast cooling scenario due to 
direct Urea process in a hyperon-mixed, pion-condensed, 
kaon-condensed, or quark-deconfined core, which gives 
softer EOS and hence a smaller maximum mass. Even 
without exotic matters, the direct Urea process may also 
become important in the core of a NS with the mass that 
is slightly less than the maximum mass. However, in the 
present paper, we exclude such relatively rare cases and 
adopt the canonical value, 1.4M Q , as the NS mass. 

On these ground, we adopt the minimal cooling sce- 
nario (Pa ge et al.l l2004). which has no enhanced cooling 
that could result from any of the direct U rea processes 
and e mploys the standard EOS, APR EOS (|Akmal et al.l 
1998). Within the minimal cooling scenario, the cool- 
ing history of a NS substantially depends on the com- 
position of the envelope, which is defined as the upper- 
most layer extending from the photosphere down to a 
boundary where the luminosity in the envelope equals 
the total surface luminosity of the star. An envelope 
contains a strong temperature gradient and has a thick- 
ness ar ound 100 meters. We adopt the cooling curves 
given in lPage et al.1 (|2004[ ) and consider the two extreme 
cases: light element and heavy element envelopes. Be- 
cause of the uncertainty in the modeling of neutron 1 So 
(i.e., spin-singiet state) Cooper paring temperature (as 
a function of neutron Fermi momentum) and proton 1 So 
paring temperature (as a function of proton Fermi mo- 
mentum), the predicted Lx distributes in a 'band' for 
each chemical composition of the envelope. When a NS 
is younger than 10 4 ' 4 yr, a light element envelope (con- 
taminated by e.g., H, He, C, or O with masses exceeding 
1O _6 M0) has a higher surface temperature, and hence 
a greater luminosity, Lx, compared to a heavy element 
envelope (contaminated by light elements with masses 
less than 10 _16 M Q ), because the heat transport becomes 
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Fig. 2. — Analytically solved gap thickness. The (blue) dot- 
ted and (red) dashed curves represent /i m = /i m (t) for a light 
and a heavy element envelopes, respectively. Uncertainties due 
to nucleon Cooper paring models, are represented by the blue and 
red shaded 'bands'. Magnetic inclination angle is assumed to be 
a = 60°. 

more efficient in the former. As the NS ages, a star with 
a light element envelope quickly loses its internal energy 
via photon emission; as a result, after 10 46 yr, it become 
less luminous than that with a heavy element envelope. 
A realistic cooling curve will be located between these 
two extreme cases, depending on the actual composition 
of the NS envelope. 

We present the solved h m for a light and a heavy el- 
ement envelope in figure [21 adopting fi^o = 3.2, which 
gives the magnetic field strength of 4.1 x 10 12 G at the 
magnetic pole, where r* = 11.6 km is used. It follows 
that the gap becomes thinner for a light element case 
(between the thin and thick dotted curves, blue shaded 
region) than for the heavy element cases (between the 
thin and thick dashed curves, red shaded region). This 
is because the more luminous photon field of a light ele- 
ment envelope leads to a copious pair production, which 
prevents the gap to expand in the trans-field direction. 
As a result, the predicted L 7 becomes less luminous for 
a light element envelope than a heavy one. 

In figure [31 we present the analytical results of L 7 ver- 
sus L S piYi as the dotted (or dashed) curve for a light (or a 
heavy) element envelope. As the pulsar spins down, L 7 
evolves leftwards. In this section, for analytical purpose, 
we are not interested in the dependence on the observer's 
viewing angle, £. Thus, in equation (QJ, we simply put 
fa = 1. It is interesting to note that L 1 little evolves 
for a light element envelope, because of h m <C 1, as dis- 
cussed after equation (fTTj) , Precisely speaking, we obtain 
h m cx i 62 and hence L 1 cx L S p in 07 (or h m cx t - 50 and 
hence L 1 cx £ sp in ' 25 ) for a light (or a heavy) element 
envelope. That is, although h m is smaller, h m increases 
more rapidly in a light element case than in a heavy 
element case, which enables a constant L 7 for a light el- 
ement envelope. At later stage, t > 10 4 years, both h m 
and L 7 increase with decreasing L sp i n for a light element 
envelope, because Lx/(kT), and hence the pair produc- 
tion rate, rapidly decreases with increasing t. On the 
other hand, for a heavy element envelope, their greater 
h m results in a monotonically decreasing L~ i with de- 
creasing i S p in . 

3. NUMERICAL EXAMINATION OF OUTER-GAP 
ELECTRODYNAMICS 

Let us develop the analytical examination and look 
deeper into a self-consistent solution by a numerical 
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Fig. 3. — Analytically solved gap luminosity with uncertainties 
due to nucleon Cooper paring models. The thick and thin dotted 
(or dashed) curves represent the evolution of the gap luminosity 
for a light (or a heavy) element envelope. Two cases of magnetic 
inclination angle, a = 45° and 60°, are depicted as labeled. The 
green filled circles designates the normal gamma-ray pulsars, while 
the blue filled squares do those detected by the gamma-ray blind 
search technique. 

metho d. To this end, we adopt the modern outer-gap 
model (|Hirotanill2011aI ) and solve the set of Maxwell and 
Boltzmann equations self-consistently and compute Eu , 
distribution functions of e s, and the photon specific 
intensity at each point in the three-dimensional pulsar 
magnetosphere. We consider not only the whole-surface, 
cooling NS emission but also the heated polar-cap emis- 
sion as the photon source of photon-photon pair pro- 
duction in the numerical analysis. The former emission 
component is given as a function of the pulsar age from 
the minimum cooling scenario, in the same manner as 
in the analytical examination, while the latter emission 
component is solved consistently with the energy flux of 
the e~'s falling on to the pulsar polar-cap surface. The 
method described below was also used in recent theo- 
retical computation o f the Crab pulsar's 7-ray emissions 
(| Aleksic et al.11201 lallbh . 

Let us present the basic equations that describe a 
pulsar outer-magnetospheric accelerator, extending the 
method first proposed for black-hole m agnetospheres 
(|Beskinlll99llHirotani fc Okamoto1 IT998). Around a ro- 
tating NS, the bac kground geometry is des cribed by the 
space-time metric ([Lense fc Thi rring 1918h 

ds 2 = g tt dt 2 +2gt v dtd(p+g rr dr 2 +gegd9 2 +g ipip dip 2 , (30) 
where 

gu = (l - ^) c 2 , gt v ee ac^ sin 2 9, (31) 

g rr = - (l- — \ , gee = -r 2 , g w = -r 2 sin 2 9; 

(32) 

M denotes the NS mass, r g ee 2GM/c 2 the Schwarzschild 
radius, and a = IQ/(Mc) the stellar angular momentum. 
At radial coordinate r, the inertial frame is dragged at 
angular frequency uj ee —gt v /g vv — 0.15f2i45r6 -3 , where 
^45 ee J/10 45 erg cm 2 , and r$ ee r*/10km. 

First, let us derive the Poisson equation for the electro- 
static potential under the space-time geometry described 
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just above. Let us consider the Gauss's law, 



An 

'- —P' 
c z 

(33) 

where V denotes the covariant derivative, the Greek in- 
dices run over t, r, 0, 



'; V~9 = V9rrge9Pi = cr 2 sin9 
and p\ = g 2 ^ — gttg^, P the real charge density. The 
electromagnetic fie lds observed by a distant static ob- 
server are given by (|Camenzindll986aT lbl) E r = F rt , Eg = 
F et , E v = F vU B r = {g tt + g tv Q)Fg v /^, B 6 = 
(gtt+gtiptyFpr/ '\/ = <?, B v = -plFre/y^g, where F^ u = 
0^A U — d„A^ and denotes the vector potential. 

We assume that the electromagnetic fields are un- 
changed in the co-rotating frame. In this case, it is conve- 
nient to introduce the non-corotational potential ^ that 
satisfies 



F, 



nF„ 



(34) 



where p = t,r, 9, <p. If Fa± + $IFa v — holds for A = r 
and 9, the magnetic field rigidly rotates with angular 
frequency fl. Imposing that the NS surface a perfect 
conductor, Fg t + QFg^ = 0, we find that the NS sur- 
face becomes equi-potential, Og^ = + HO^ = 0. 
However, in a particle acceleration region, the magnetic 
field does not rigidly rotate, because FAt + QFa v devi- 
ates from 0. The deviation is expressed in terms of "J, 
which gives the strength of the acceleration electric field, 



" B B y 



m v ) = § • (-V*), (35) 



which is measured by a distant static observer, where the 
Latin index i runs over spatial coordinates r, 0, ip. 

Substituting equation (f34|) into (|33l) . we obtain the 
Poisson equation for 



--d„ 



g^g vv d u ^)=A7r{p-p G j), (36) 



where 



Pgj 



(37) 



denotes the general relativistic Goldreich- Julian charge 
density. If p = pgj holds everywhere, E\\ vanishes in 
the entire region, provided that the boundaries are equi- 
potential. However, if p deviates from pgj m & n y region, 
2£|| appears around the region with differentially rotat- 
ing magnetic field lines. In the higher altitudes, r ^> r g , 
equa tion (I3T1) reduces to the special-relativi stic expres- 
sion (|Goldreich. fc JulianlH969l : IMestell[l97l , 



ft B (flxr)'(Vx B) 



2nc 



Aire 



(38) 



which is commonly used. 

Instead of (r,8,ip), we adopt the magnetic coordinates 
(s,0*, </?*), where s denotes the distance along a mag- 
netic field line, 8* and represents the magnetic co- 
latitude and the magnetic azimuthal angle, respectively, 
of the point where the field line intersects the NS sur- 
face. Defining that 9^,-0 corresponds to the magnetic 
axis and that if* — to the latitudinal plane on which 



both the rotation and the magnetic axes reside, we ob- 
tain the Poisson equation {Hirotani 2006a]13), 



Pi 



[g^dl+g^'dl+g^d, 



+2g st) '0 s 0g, + 2g°**»d e ,0 tp , + 2g^ s Vt s ) * 
-{A s s +A e '0g, +A*-d v ,)y = 47r(p-p G j), (39) 
where 

,0x l 0x j ^„0x l dxi ggOx 1 0x j 



g tJ = 9» v 



Ox^ dx v ' Or Or 
k Ox 1 ' dx^' 
Pi dip ' 
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(40) 




the coordinate variables are x 1 = r, x 2 = 9, x = tp, 
x 1 = s, x' 2 = 9*, and x 3 = </?*. Note that this formalism 
is applicable to arbitrary magnetic field configurations, 
and also that equation ([35]) gives E\\ — — (d^ /Os)e,, Vt , . 
Causality requires that a plasmas can co-rotate with 
the magnetic field only in the region that satisfies kg = 
gtt + 2g tip fl + g w & 2 > (|ZnaielJH977l : iTakahasi et all 
1990). That is, the concept of the light cylinder is gen- 
eralized into the light surface on whi ch fcn vanishes. The 
effect s of magnetic field expansi on (Sch arlemann et all 
119781 : iMuslimov fc TsvganlTl99l are contained in the 
coefficients of the trans-field derivatives, g e * e * , g 6r(p * , 
g<p,tp*_ T n wna t foll ows, we adopt the vacuum, rotat- 
ing dipole solution ([Cheng et al.l |2000T) to describe the 
magnetic field configuration. 

Second, let us consider the particle Boltzmann equa- 
tions. At time t, position r, and momentum p 7 they 
become, 



^± + v ■ VN ± + (qE+"xB 



0N ± 
Op 



= S±(t,r,p), 
(42) 

where N+ (or 7V_) denotes the positronic (or electronic) 
distribution function; v =p/(m c 7), m c refers to the rest 
mass of the electron, and q the charge on the particle. 
The Lorentz factor is given by 7 



1/y/l - (\v\/c) 2 . The 
collision term S+ (or S- ) consists of the terms that rep- 
resent the appearing and disappearing rates of positrons 
(or electrons) at r and p per unit time per unit phase- 
space volume. Since we are dealing with high-energy 
phenomena, we consider wave frequencies that are much 
greater than the plasma frequency and neglect the col- 
lective effects. 

It is noteworthy that the particle flux per magnetic 
flux tube is conserved along the flow line if there is 
no particle creation or annihilation. Thus, it is con- 
venient to normalize the particle distribution functions 
by the Goldreich- Julian number density such that n± = 
(N±)/(SlB/2Trce), where () denotes that the quantity is 
averaged in a gyration. Imposing a stationary condition 
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d / dt+Vtd / dip = in the co- rotating frame , we can reduce 
the p article Boltzmann equations into (Hir otani et al.l 
2003) 



ccosx- 



dru. 



Os 



P- 



,dn ± ,dn± 
"5 r — — OiC + Op, 

Op o\ 



(43) 



where the upper and lower signs correspond to the 
positrons (with charge q = +e) and electrons (q = — e), 
respectively, x denotes the pitch angle of gyrating parti- 
cles, and p = \p\ = m c cy/ 7 2 — 1. Since pair annihilation 
is negligible in a pulsar magnetosphere, S± consists of the 
IC scattering term, Sic, and the pair creation term, S p , 
in the right-hand side. When a particle emits a photon 
via synchro-curvature process, the energy loss (~ GeV) 
is small compared to the particle energy (~ 10 TeV); 
thus, it is convenient to include the back reaction of the 
synchro-cu rvature emission as a friction term in the left- 
hand side (Hirotani 2006a). In this case, the character- 
istics of equation (|4"3")l in the phase space are given by 



p = qE\\ cosx 



Psc 



X 



qEn sinx d(\nB^ 2 ) 
he smx, 

V OS 



(44) 
(45) 



where the synchro-cu r vature radiation force , Psc/c 
(ICheng fc Zhand H99l IZhang fc Chengl li~99l . is in- 
cluded as the friction; the particle position s is related 
with time t by s = ds/dt = ccosx- For outward- (or 
inward-) migrating particles, cosx > (or cosx < 0). If 
En =0, particles will be reflected by magnetic mirrors, 
which are expressed by the second term in the right-hand 
side of equation (|45l If we integrate n± over p and Xi an d 
multiply the local GJ number density, QB/2irce, we ob- 
tain the spatial number density of particles. Therefore, 
we can express the real charge density p as 



VLB 

2^c 



(n+ - n„)djdx + ph 



(46) 



where n± are a function of s, 0», <p*, 7, x; Pion refers to 
the charge density of ions, which can be drawn from the 
NS surface as a space-charg e-limited flow (SCLF) by a 
positive E\\ (|Hirotanill2006a| ). 
In equation (j43|) . the collision terms are expressed as 



-l<7 



'IC : 



^c(( e «> 



■If 



n±djdx, 



+ '52vic(li,l,Vc) // n±d-ydx 



(47) 



and 



lixce f [In 
S P = -jTg / ("7-8 + "77) / j—dtidv, (48) 



where (e„) = (e n _i + e n )/2 represents the typical pho- 
ton energy within the interval [m e c e„_i,m e c e n ], p, c the 
cosine of the collision angle between the particles and 
the soft photons, I v the specific intensity of the radia- 
tion field, Q the solid angle into which the photons are 
propagating. To compute S p at each point, we have to 
integrate /„/ (hv) in all directions to calculate the differ- 
ential photon number flux, J I v / (hv)duj. The absorption 



coefficients for magnetic and photon-photon pair produc- 
tion processes, are given by 



e 2 /He m e 



-B'_ 



exp 



1 



3e 7J B' 



4.4 h 

X(5(7 - 7o)<5(x - Xo), 
l-cos(9 r f°° , dF s d 2 a. 



(49) 
(50) 



where B' 



Bsm6 c /B c 



de s djdx 
explicit expression of 



77n and d 2 a^y / d^dy are given in the literature (e.g 
Hirot ani et all (|2003l )) 



In 



equation (|4"9")l . B'± con- 
tains the collision angle, 9 C , between the photon and 
the magnetic field. In equation (j5"0")l . 9 C does that be- 
tween the two photons. The IC redistribution function 
7?7c( e 7) 7' Mc) represents the probability that a particle 
with Lorentz factor 7 up-scatters photons into energies 



between 



and 



"(e 7 + de-y) per unit time when 



the collision angle is cos /i c , where e 7 = hv/{m e c ) 
refers to the dimensionless photon energy. On the other 
hand, f7ic(7i>7/> A*c) represents the probability for a par- 
ticle to change its Lorentz factor from 7, to 7/ in a sin- 
gle scattering. We thus obtain ?7ic (7i j 7/ 7 /^c ) = 7 ?7c(7« — 
7/,7i, He) by energy conservation. 

Third, let us consider the radiative transfer equation. 
The variation of specific intensity, I„, along a ray is de- 
scribed by the radiative transfer equation, 

^jj- = -a u I v + j u , (51) 

where I refers to the distance along the ray, a v and j v the 
absorption and emission coefficients, respectively. Both 
a v and j v are a function of I, photon energy E-y, and 
propagation direction {k e ,k v ), where k* 1 denotes the pho- 
ton momentum four vector (ji = t, r, 9, if). 

Equation (|5l"|) can be solved if we specify the photon 
propagation in the curved space time. The evolution of 
momentum and position of a photon is described by the 
Hamilton- Jacobi equations, 

dk L __dh ^ ^ 



(53) 





dh 


dkg 


dk t 


° dl 


dr ' 


c ~dT 


de ' 


dr 


dh 


dO 


_ dh 


c di = 


dk r ' 




dkg' 



Since the metric (eqs. [30] 



j) is stationary and ax- 
isymmetric, the photon energy at infinity k t and the az- 
imuthal wave number —k v are conserved along the ray. 
When a particle is rotating with angular velocity <p and 
emit a photon with energy -Eqocab h and —k v are re- 
lated to these quantities by the redshift relation, E\ oca \ = 
[dt / dr)(kt + half), where dt/dr is solved from the defini- 
tion of the proper time, (dt / dr) 2 (gtt + 2gt v ip + g vv ^> 2 ) — 
1. The dispersion relation fc M fc M = 0, which is quadratic 
to k^s (/1 = t,r,9,ip), gives Hamiltonian k t in terms of 
r, 6, k r , kg, and k v . Thus, we have to solve the set of 
four ordinary differential equations (|52|) and (f53|) for k r , 
kg, r, and 6. When photons are emitted, they are highly 
beamed along the particle's motion; thus, the initial con- 
ditions of (k r ,k ,k v ) are given by the instantaneous par- 
ticle's velocity measured by a distant static observer. It 
may be helpful to give the initial conditions of ray trac- 
ing in the limit r ^> r g and 7 > 1, because photons are 
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emitted mostly in the outer magnetosphere. In this limit, 
the orthonormal components o f the e + ' s and e~'s instan- 
taneou s velocity are given by (|MesteHll985t iCamenzindl 
H986allbT l 

v r B r v § . B 9 , w 

— = /v — , — = /v— , — = /v + (54) 

C C C C C C ^LC 

in the polar coordinates, where 




B 



p = y (B r ) 2 + (B 9 ) 2 . The upper (or the lower) sign of 
/ v (eq. [33]) corresponds to the outward (or the inward) 
particle velocity. Note that B 1 ^ < holds in ordinary 
situation. 

Fourth and finally, let us impose appropriate bound- 
ary conditions to solve the set of Maxwell (i.e., Poisson) 
and Boltzmann equations. We start with considering the 
boundary conditions for the elliptic type equation (|36[) . 
We define that the inner boundary, s = 0, coincides with 
the NS surface, on which we put ^ = for convenience. 
The outer boundary, s = s ou t (0» ,<p*), is defined as the 
place where Eu changes sign near the light cylinder. 
Its location is solved self-consistently as a free-boundary 
problem and appears near the place where d{pa / B) / ds 
vanishes due to the flaring up of the field lin es towards 
the rotation axis (eq. [68] of IHirotani (2006a)). At each 
<p», the lower boundary 8* = 0™ ax (^>*) is assumed to co- 
incide with the last open field line, which is defined by 
the condition that y/—g rr B r sin 8 + ^J—geeB 9 cos 8 = 
is satisfied at the light cylinder. We can compute the 
potential drop along each field line, A4 f (0», (£>*) = \&(s = 
0) - *(s = s out ) = = s out ), by integrating E\\ 

along the field line. The maximum value of Aif/(9*,(p*) 
is referred to as A9 mBX . 

Let us consider the upper boundary, 0J" n = 
0™ ln (s, 93*). At each Eu vanishes on the last- 

open field line, 6>* = 0™ ax (i. e., at the lower bound- 
ary), and increases with decreasing 0* in the latitudi- 
nal direction (i. e., towards the magnetic axis). The 
acceleration field Eu peaks around the middle-latitudes, 

0» » [9? ax + 8f n }/2 = (1 - h m /2)8™ a *, and turns to de- 
crease towards the upper boundary. At some co-latitude, 
En eventually decreases below 0.01A>I' max /'Ci7Lc; we de- 
fine this co-latitude as the gap upper boundary, 0™ m = 
0™ ln (s, y>»). To specify the magnetic field line at each 
(p», it is convenient to introduce the dimensionless co- 
latitude, 

h=l- 9*/9™ ax (ip*). (56) 

For example, h — specifies the last-open field line, and 
ft, = 0.1 does the field line having the foot point on the 
polar-cap (PC) surface at 0* = O.90" iax at each magnetic 
azimuthal angle <p*. At the upper boundary, we obtain 



h = Ms, = °f n (s, <p*)/9? ax (<p*). 



(57) 



To solve the Poisson equation, we put $ = at h = h m , 
or equivalently, at 9* = 8™ m (s,tp*). If h m <C 1, the gap 
becomes thin, whereas h m ~ 1 indicates that the gap 
is threaded by most of the open-field lines. Note that 



the open field lines cross the PC surface at magnetic co- 
latitudes 0™ = (1 - h m )9f ax < 8* = (1 - /i)0„ max < 
0™ ax (i.e., h m > h > 0), and that 9* = (i.e., h = 1) 
corresponds to the magnetic axis. 

We also have to consider the boundary conditions for 
the hyperbolic type equations (|43|) and ([ST]) . (Eq. [51] 
itself is an ordinary differential equation; however, it is 
equivalent to solving the Boltzmann equation of photon 
distribution function.) We assume that e^s and photons 
are not injected into the gap across either the inner or 
the outer boundaries. However, if the created current 
becomes greater than the GJ value, a positive Eh arises 
at the NS surface to draw ions from the surface as a 
SCLF until Eh almost vanishes at the surface. 

To sum up, we solve the set of partial and ordinary dif- 
ferential equations ([36]) , ([43]) , and (|5ip under the bound- 
ary conditions described in the three foregoing para- 
graphs. By this method, we can solve the acceleration 
electric field Eu, particle distribution functions n±, and 
the photon specific intensity I v (from hv = 0.005 eV to 
50 TeV) , at each position in the three-dimensional mag- 
netosphere of arbitrary rotation-powered pulsars, if we 
specify P, [i, a, and kT, where the surface tempera- 
ture kT is necessary to compute the photon-photon pair 
production through the differential flux dF s /de s in equa- 
tion (|50p . We adopt the minimum cooling scenario in the 
same manner as in § [5] 

In figure [4] we plot the result of L 7 as a function of 
L sp in as the dash-dotted (or solid) curve for a light (or a 
heavy) element envelope, where /X30 = 3.2 is adopted in 
the same manner as in the analytical examination. It fol- 
lows that these numerical solutions are consistent with 
the analytical ones, and that L 1 decreases slowly until 
10 4 5 years. The physical reason why L 1 increases with 
decreasing i sp i n at t > 10 4 years for a light element enve- 
lope, is the same as described at the end of § El A realistic 
NS will have an envelope composition between the two 
extreme cases, light and heavy elements. Thus, the ac- 
tual L 7 's will distribute between the red solid (or dashed) 
and the blue dash-dotted (or dotted) curves. However, 
after approaches L Rn \ n (thin dashed straight line; see 
iWang fc Hirotanil ([20111 ) for the death line argument), 
the outer gap survives only along the limited magnetic 
field lines in the trailing side of the rotating magneto- 
sphere because of a less efficient pair production; as a 
result, L.y rapidly decreases with decreasing L spin . For a 
smaller a, even for a light element envelope, L 7 mono- 
tonically decreases as the dash-dot-dot-dot curve shows, 
because the gap is located in the higher altitudes, and 
because the less efficient pair production there prevents 
the produced electric current to increase with decreasing 
age around t ~ 10 4 5 years. 



4. EXPONENTIAL CUTOFF ENERGY 

It is also worth examining the cutoff energy, Ecutoff , of 
the gamma-ray spectrum. We plot -E C utoff &s a function 
of the magnetic-field strength at the light cylinder, -Blc, 
in figure O In the analytical computation (dash, dotted, 
and dash-dot-dot-dot curves), only the curvature pro- 
cess ([Hirotanil 1201 lal ) from the particles produced and 
accelerated in the outer gap, is considered as the pho- 
ton emission process, and no photon absorption is con- 
sidered. In the numerical computation (solid and dash- 
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Fig. 4. — Same figure as figure [3] but numerical results are 
plotted as dash-dotted and solid curves for a light and a heavy 
element envelope, respectively. For comparison, a numerical solu- 
tion for a = 45° and a light element envelope, is also plotted as 
the dash-dot-dot-dot curve. To avoid complications, we adopt the 
geometrical mean of the upper-bound Lx and lower-bound Lx of 
cooling curves for each chemical composition, instead of depicting 
a 'band', which reflects the uncertainties in the nucleon Cooper 
pairing models. 

dotted curves), on the other hand, synchro-curvature and 
inverse-Compton processes from the particles not only 
produced and accelerated in the gap but also cascaded 
outside the gap, are considered as the emission processes, 
and both the photon-photon and magnetic pair pro- 
duction processes are taken into account when comput- 
ing the photon propagation. Therefore, for very young 
pulsars, a strong photon-photon absorption (and the 
subsequent production of lower-energy photons via syn- 
chrotron and synchrotron-self-Compton processes) takes 
place even in the higher altitudes. As a result, spectrum 
can no longer be fitted by power-law with exponential- 
cutoff functional form. Thus, for t < 10 4 years, the fit- 
ted cutoff energies are not plotted for numerical solutions 
(i.e., solid and dash-dotted curves). 

It follows that the present outer gap model explains the 
observed tendency that -©cutoff increases with increasing 
Slcj because the Goldreich- Julian charge density (in the 
gap) increases with increasing L?lc- It also follows that 
©cutoff is regulated below several GeV, because copious 
pair production leads to h m <C 1 for young pulsars (i.e., 
for strong -Blc)- 

5. FLUX CORRECTION FACTOR 

Finally, let us investigate the flux correction factor, /q. 
To infer the 7-ray luminosity, L 1 = Anf^F 1 d 2 , from the 
observed flux, Ky, one has conventionally assumed that 
the flux conversion factor is approximately unity, /n ~ 1, 
where d denotes the distance to the pulsar. However, 
now we can compute /n explicitly as a function of the 
observer's viewing angle, £, using the three-dimensional 
numerical solution. In figure [HI /n(C) is depicted for 
heavy and light element envelopes at 10 4 and 10 5 years. 
It follows that the error of /n = 1 is kept within a fac- 
tor of 3 with a probability greater than 50%. When the 
pulsar is young, Eu is highly screened in the middle and 
lower altitudes; as a result, most 7-rays are emitted from 
the higher altitudes to appear within the observer's view- 
ing angle 46° < C < 53° (i.e., 0.6 < cosC < 0.7). How- 
ever, as the pulsar ages, more 7-rays are emitted from 
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Fig. 5. — Cutoff energy of 7-ray spectrum. The dashed and 
dotted curves denote the cutoff energies obtained analytically for 
heavy and light element envelopes, respectively, whereas the solid 
and dash-dotted curves do those obtained numerically for heavy 
and light element envelopes, respectively. For these four cases, 
a = 60° is assumed. For comparison, an analytical solution for 
a = 45° and a light element envelope, is also plotted as the dash- 
dot-dot-dot curve. 
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Fig. 6. — Flux correction factor, /q, as a function of the ob- 
server's viewing angle, £, with respect to the rotation axis. The 
(red) thick solid and dashed curves denote /n at pulsar age I0 4 and 
10 5 years, respectively, for a heavy element envelope, whereas the 
(blue) thin solid and dashed curves do those at 10 and 10 5 years, 
respectively, for a light element envelope. The abscissa denotes the 
cosine of £; thus, C distributes randomly along the abscissa with 
uniform probability. For all the cases, a = 60° is assumed. 

the middle and lower altitudes, resulting in a stronger 
flux near the rotational equator, 66° < ( < 73°. 

6. DISCUSSION 

To sum up, a light element envelope approximately 
corresponds to the lower bound of the (observationally 
inferred) gamma-ray luminosity of rotation-powered pul- 
sars, whereas a heavy element one to the upper bound. 
The scatter of the intrinsic gamma-ray luminosity is 
physically determined by the magnetic inclination an- 
gle, a, and the envelope composition. The cutoff energy 
of the primary curvature emission is kept below several 
GeV even for young pulsars, because the gap trans-field 
thickness, and hence the acceleration electric field, is sup- 
pressed by the polarization of the produced pairs in the 
lower altitudes. 

To convert the observed 7-ray flux into luminosity, Ly, 
one has conventionally assumed fn = 1. For example, 
the error bars of the observational data points in figurelH 
do not contain any uncertainties incurred by /q. Never- 
theless, if a and £ can be constrained, we can estimate 
Ly more accurately, by applying the present quantita- 
tive outer-gap calculations. It is noteworthy that L 7 's 
given in figures [3] and 0] little depend on the NS magnetic 
moment, fi. This is particularly true for a light element 
case, which has h m <C 1, by the reason described after 
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equation (1111) . What is more, with an additional determi- 
nation of d (e.g., by parallax observations), we can infer 
the composition of individual NS envelopes, by using the 
constrained flux correction factor, /n (fig. EJ. We hope 
to address such a question as the determination of a and 
C, and hence /o, for individual pulsars, by making an 'at- 
las' of the pulse profiles and phase-resolved spectra that 
are solved from the basic equations in a wide parameter 
space of P, fi, T, a, and £, and by comparing the atlas 
with the observations. 
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